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Abstract 

We simulate the transition from amphfied spontaneous emission (ASE) to lasing in random 
systems with varying degrees of mode overlap. This is accomplished by solving the stochastic 
Maxwell-Bloch equations with the finite-difference time-domain method. Below lasing threshold, 
the continuous emission spectra are narrowed by frequency-dependent amplification. Our simula- 
tions reproduce the stochastic emission spikes in the spectra. Well-defined peaks, corresponding to 
the system resonances, emerge at higher pumping and are narrowed by stimulated emission before 
lasing takes place. Noise tends to distribute pump energy over many modes, resulting in multi- 
mode operation. Well above the lasing threshold, the effects of noise lessen and results become 
similar to those without noise. By comparing systems of different scattering strength, we find that 
weaker scattering extends the transition region from ASE to lasing, where the effects of noise are 
most significant. 
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I. INTRODUCTION 



A random laser has two elements: an active material that spontaneously emits light and 
amplifies it via stimulated emission, and a disordered medium that partially traps light by 
multiple scattering. Lasing occurs when the loss due to absorption and leakage of light 
through open boundaries is compensated by light emission and amplification inside the 
medium. Stronger scattering increases the trap time of light and lowers the leakage loss. 
Letokhov discussed the diffusion process with gain that can lead to lasing with non-resonant 
feedback in the 1960s Early experiments, e.g., on dye solutions containing microparticles 
2! , showed a dramatic narrowing of the emission spectrum and a rapid increase of emission 
intensity at the frequencies around the maximal gain where the threshold condition is met. 
In contrast to the smooth and relatively broad lasing spectra, later experiments illustrated 
multiple sharp peaks of laser emission from semiconductor powder and disordered polymers 
[3I, Q . Those spectral peaks result from interference of scattered light in the random media. 
Although interference is not required for lasing action, it reduces light leakage at certain 
frequencies [5| and facilitates lasing by lowering the threshold (gain = loss). Thus lasing 
occurs at those frequencies, producing emission with high first-order coherence (narrow 
spectral width) and second-order coherence (suppression of photon number fiuctuations 



in single modes) [6]. In addition to the reproducible lasing peaks, stochastic spikes were 
observed in single-shot emission spectra with pulsed excitation [l-lO]. These spikes are 
completely different in frequency from shot to shot, and are attributed to strong amplification 
of spontaneously emitted photons along long paths. 

The rich phenomena of random lasers have not been well understood so far. The diffusion 
model including gain can describe the narrowing of broad emission spectra but not the 
appearance of discrete lasing peaks since it neglects the interference effects. Semiclassical 
laser theory, based on the Maxwell's equations, can predict the lasing peaks but not 

the stochastic emission spikes because it does not take into account the spontaneous emis- 
sion. The spectral width of lasing modes cannot be calculated either without spontaneous 
emission. We do not know how the laser linewidth compares to the frequency spacing of 
lasing modes. If the former is larger than the latter, the lasing peaks are indistinguishable 
no matter how fine the spectral resolution is. A diffusive or ballistic system, especially 
with higher dimensionality, contains a huge number of resonances that overlap spatially 
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and spectrally. Although these resonances have similar lasing thresholds, they cannot all 



lase simultaneously because of gain depletion Jl3| 



have been observed experimentally 



0, y, c 



A large fluctuation of lasing spectra 



2l|. The lasing modes are sensitive to small 



perturbations and noise, e.g., fluctuation of pump pulse energy, spatial variation of pump 
intensity, etc. [l4|- In addition to the extrinsic noise, the number of spontaneously emit- 
ted photons participating in the buildup of laser emission in any mode may fluctuate from 



shot to shot, leading to variations of lasing peak height [20|, |22|. Such intrinsic fluctuations 
are missed by semiclassical laser theory. Since random laser thresholds are usually higher 
than conventional laser thresholds due to weaker optical conflnement, stronger pumping is 
required, making the amplifled spontaneous emission (ASE) stronger. Semiclassical laser 
theory, which neglects ASE, cannot capture the transition from the ampliflcation of spon- 
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taneous emission to lasing oscillation that has been observed experimentally 

Note that spontaneous emission only contributes to part of the intrinsic noise, which 
also includes the fluctuations induced by atomic dephasing, pumping, optical leakage, etc. 
Since intrinsic noise plays an essential role in random lasing behavior, it must be treated 
properly. There have been signiflcant advances in theoretical studies on photon statistics of 



random lasers and ampliflers 



24j428l|. Most of them are based on full quantum treatments 



of noise in the modal description. For a random system, the mode structures are complex 
and unknown a priori. Thus it is desirable to introduce noise without prior knowledge of 
modes. Some previous studies based on light diffusion and random walks 29| do not need 
mode information, but they ignore light interference that is essential to the formation of 
resonant lasing modes. 

In this paper we incorporate intrinsic noise into the numerical simulation of random lasers. 
The numerical method is based on the flnite-difference time-domain (FDTD) formulation 
we recently developed to study the effects of noise on light-atom interaction in complex 
systems without prior knowledge of resonances |30|, ImI]- The interference effects and the 
openness of the system are fully accounted for with Maxwell's equations and absorbing 



boundary conditions 



32j | . The incorporation of the Bloch equations for the density of states 



of atoms simulate the dynamics of atoms and their interaction with light 



33| . In the many- 



atom and many-photon limit the quantum fluctuations can be simulated by classical noise 



terms 



34| . Based on the fluctuation-dissipation theorem, we consider noise associated with 



three dissipation mechanisms for atoms (described in detail in j3l|) (i) dephasing events. 
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(ii) excited state decay, (iii) incoherent pumping (from ground state to excited state). Noise 
related to field decay is negligible because the photon energy at visible frequencies is much 
larger than the thermal energy at room temperature. At higher temperatures or longer 
wavelengths, where this noise becomes significant, it can be incorporated into the FDTD 
algorithm following the approach we developed in 30 1. 

Here we study the effects of intrinsic noise on the steady-state properties of random 
lasers in one dimension (ID). Results from systems of different scattering strengths are 
presented which probe varying degrees of light leakiness and spectral mode overlap. For the 
first time, we are able to simulate the transition from ASE to lasing using the stochastic 
Maxwell-Bloch equations. Stochastic emission spikes are reproduced with similar statistics 

nn 

to the experimental data reported previously [9|, 110[. The spectral width of the broad ASE 
peak is calculated as a function of pumping rate. It displays a dramatic decrease with 
increasing pump level as seen experimentally . The linewidths of individual lasing modes 
are also computed and compared to the Schawlow-Townes linewidth of single mode lasers. 
A comparison of the results of simulations with noise to the simulations of the same active 
system without noise illustrates that noise effects are strongest in the transition regime from 
ASE to lasing. 

This paper is organized as follows. In Sec. HT] information of the random systems studied 
here is provided. An analysis of resonances in these systems without gain is carried out in 
Sec. Illli In Sec. IIVI the FDTD formulation for the stochastic Maxwell-Bloch equations are 
given and some numerical issues related to noise are discussed. Results of calculations using 
the Maxwell-Bloch equations both with and without noise are presented for random systems 
with spectrally overlapping resonances in Sec. |V] and with non-overlapping resonances in 
Sec. I VII Our main conclusions are drawn in Sec. IVIII 



II. RANDOM SYSTEMS 



Two ID random systems are considered here with different degrees of spectral overlap 
of resonances. Both consist of = 41 layers. The dielectric layers with index of refraction 
rii > 1 alternate with air gaps {n2 = 1) resulting in a spatially modulated index of refraction 
n{x). The scattering strength is varied by adjusting the index contrast An = ni/n2 — 1. 
The system is randomized by specifying different thicknesses for each of the layers as di^2 = 
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{di,2) (1 + where {di) and (^2) are the average thicknesses of the layers, < < 1 
represents the degree of randomness, and C is a random number uniformly distributed in 
(-1,1). The average thicknesses are {di) = 100 nm and (^2) = 200 nm, giving a total average 
length of (L) = 6100 nm. The grid origin x = is at the left boundary of the structure, and 
the length of the random structure L is normalized to (L). The degree of randomness is set 
to ?7 = 0.9 and the index of refraction outside the random media is no = 1. 

The degree of mode overlap is adjusted by the refractive index rii of the dielectric layers. 
The Thouless number g, which reveals the amount of spectral overlap of resonances of these 
random systems, is given by the ratio of the average resonance decay rate to the average 
frequency spacing g = {k.j) / (Ak). In the first case, ni = 1.05 {An = 0.05), g = 1.0 and the 
resonances overlap in frequency. In the second case, rii = 1.25 {An = 0.25), g = 0.5 and the 
resonances are fairly well separated. The localization length ^ is obtained from the variation 
of transmission over the system length, and averaged over the wavelength range of interest 
(400 nm to 1200 nm). (^) = 340 /im for the first case, and (^) = 13 fim for the second one. 
Since (C,) is much larger than the system length L, both systems are far from localization 
threshold. Figure [H^a) shows the transmission spectra T{k) of both systems. Resonance 
peaks are clearly narrower and better separated for the second system with g = 0.5. 



III. RESONANCES OF THE PASSIVE SYSTEM 



We calculate the resonances of the two systems in the absence of gain or absorption 
using the transfer matrix method. Because the system is open, light can escape through 
the boundaries. To satisfy the conditions that there are only outgoing waves through the 
boundaries, the wave vectors must be complex numbers, k = k + iki. The real part k 
corresponds to the mode frequency u, k = u/c, where c is the speed of light in vacuum. 
The imaginary part ki < 0; its amplitude is proportional to the decay rate 7 of the mode 
due to light leakage. The resulting field distributions associated with the solutions for these 
boundary conditions are the quasimodes of the passive system. Figure [D^b) plots the electric 
field intensity distributions of representative quasimodes in the two systems studied here. 
With g = 0.5 the spatial distribution of electric field intensity is more concentrated inside 
the system than with g = 1.0 where intensity distribution is concentrated on the boundaries 
of the system. 
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FIG. 1: (a) Transmission spectra T{h) of passive random systems with g = 0.5 (solid line) and 
5 = 1.0 (dashed line). The gain curve (dotted line) for the Maxwell-Bloch simulations is also 
shown, (b) Intensity distribution of a representative quasimode in a random system with g = 0.5 
(solid line) and 5 = 1.0 (dashed line). 

Figure |2] shows all quasimodes of the two systems in the complex-/c plane within the 
wavelength range of interest. For (7 = 1.0 the separation of decay rates between neighboring 
modes is quite small, leading to significant spectral overlap of modes. With gain included, 
the lasing thresholds of most modes should be very similar. However, the simulations of 
lasing in the following sections include a frequency- dependent gain curve centered at ka 
(shown as the vertical line in Fig. |2]). A balance of lower decay rate (smaller amplitude of 
ki) and higher gain [k closer to ka) selects the modes that are amplified the most, e.g., the 
two modes circled in Fig. |2l 

For the system with g = 0.5 in Fig. Ej the resonances have smaller amplitudes of ki and 
thus lower decay rates. This is a result of greater confinement of light due to the higher 
index contrast. Furthermore, decay rates are more separated in general. As suggested by 
the narrow peaks around k = 10.5 /xm~^ in the transmission spectrum in Fig. [H^a), the two 
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FIG. 2: Frequencies k and decay rates ki of quasimodes in two random systems with g = 1.0 (+'s) 
and g = 0.5 (x's). The vertical dashed gray hne marks the atomic transition frequency ka in 
the Maxwell-Bloch simulations. The two strongest lasing modes found in the following section are 
bounded by (circles) for g = 1.0 and (squares) for g = 0.5. The circled modes have smaller decay 
rates than the modes nearest ka- 

modes nearest ka have relatively small decay rates. They are also fairly well separated from 
neighboring modes. Thus, these two modes should have the lowest lasing thresholds. 



IV. FDTD SIMULATION OF STOCHASTIC MAXWELL-BLOCH EQUATIONS 

We consider two-level atoms uniformly distributed over the entire random system. The 
atomic transition frequency is set to ka = 10.5 fim~^, the corresponding wavelength Aq = 600 
nm. The lifetime of atoms in the excited state Ti and the dephasing time T2 are included 
in the Bloch equations for the density of states of atoms. The spectral width of the gain 



regime is given by Aka = (l/Ti + 2/T2)/c [35[. We set Ti = 1.0 ps. The values of T2 



are chosen such that the gain spectrum spans ten quasimodes of the passive system, i.e., 
Aka = 10 (Ak). The average frequency spacing (Ak) is slightly different for the two cases 
studied. For g = 0.5, T2 = 1.4 fs and Aka = 4.7 /im"^. For g = 1.0, T2 = 1.3 fs and 
Aka = 5.0 yum~^, as shown by the dotted line in Fig. [U^a). We also include incoherent 
pumping of atoms from level 1 to level 2. The rate of atoms being pumped is proportional 
to the population of atoms in level 1 pu, and the proportional coefficient Pr is called the 
pumping rate. 
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To introduce noise to the Bloch equations, we used the stochastic c-number equations 



that are derived from the quantum Langevin equations in the many-atom hmit |3J]. The 
noise sources in these equations are from both the dissipation of the system and the nonhn- 
earity in the Hamiltonian. The latter represents the nonclassical component of noise, giving 
rise to nonclassical statistical behavior. Since we are interested in the classical behavior 
of macroscopic systems, such as ASE and lasing, we neglect the nonclassical noise in our 
simulation. The classical noise results from the decay, dephasing and pumping of atoms, as 
dictated by the fluctuation-dissipation theorem. The amplitude of classical noise accompa- 
nying the field decay is proportional to the square root of the thermal photon number. At 
room temperature the number of thermal photons at visible frequencies is negligible, thus 
the noise related to the field decay is ignored here. 

The stochastic simulations solve for the atomic population of the excited states P22 and 
the atomic polarization pi = pi2 + P21 and p2 = i{pi2 — P21) which couple to Maxwell's 
equations. The stochastic equations solved at each grid point in space are 

^^^^^'^^ =ckap2ix, t) - ^pl{x, t) + Tiix, t) 

at T2 

= - cKp,{x,t) + ^E,{x,t) (2p22(x,t) - N,) 
at n 

- —P2{x,t)+T2{x,t) 
^2 

dp22{x,t) . ,^ . ,^ 1 . ,^ 
= - -^E,{X,t)p2[x,t) - —p22{X,t) 

+ ^(Ar,-p22(x,t))+r22(x,t), (1) 
J-l 

where is the electric field, 7 is the dipole coupling term, is the number of atoms per 
grid cell, and the noise terms 

Ti{x,t) =2^i{t)^JpP22{x,t) 
T2{x,t) = - 2^2{t)^-fpP22{x,t) 
r22{x,t) =Ut)VP22{x,t)/T, + Pr.pii{x,t)/Ti, (2) 

where 7^ = I/T2 — l/2Ti. The terms are real, Gaussian, random variables with zero mean 
and the following correlation relation 

{Ut)Ut')) = SjkS{t-t'), (3) 
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where j,k = 1,2,3. We assume T2 ^ Ti, and the pump fluctuations in Fi and r2 are 
neglected because they are orders of magnitude smaller than the noise due to dephasing. 
The resulting Maxwell-Bloch (MB) equations are solved through a parallel FDTD imple- 
mentation with the spatial grid step Ax = 1.0 nm and the temporal step At = 1.7 x 10^^® 
s. In cases where noise is not included, the system is excited by a Gaussian-sinusoidal pulse 
of center frequency ko = ka and spectral width Ako = Aka- 

An issue concerning the simulation with noise arises when pn or P22 is close to zero. 
To keep the atomic populations in both levels positive for a large range of pumping rates, 
we set the system initially at the transparency point, i.e., p^it = 0) = P22 — Pn = 0. 
Moreover, we assume the atomic density is large, Natom/V = 4.3 x 10"'^'^ cm~^ [36] • Small 
variations to the initial population do not affect the final steady-state results. Furthermore, 
the high frequency components of the noise excite the modes resonating within single air gaps 
sandwiched between dielectric layers of index ni > 1. These high frequency contributions 
are ignored completely by considering only the electromagnetic fields within the wavelength 
range 400 nm < A < 1200 nm. 

With noise terms included in the Maxwell-Bloch equations, all quantities fluctuate in 
time. Eventually their values averaged over small time windows are nearly constant. By 
comparing the spectra of output light taken over different temporal ranges up to t = 267 ps, 
we find a steady state is reached by 16.6 ps for all pumping rates considered here. Hence, 
the output spectra obtained after 16.6 ps represent the steady state behavior. The output 
field is sampled at the grid point x = L at the right boundary of the random system. The 
results from this point are identical in character to those from any point outside the system 
and before the absorbing boundary layer. 

V. ASE AND LASING IN A SYSTEM WITH OVERLAPPING RESONANCES 
A. Input-output relation 

Starting from the random system with g = 1.0, we investigate the transition from sponta- 
neous emission to ASE and to laser emission by examining the dependence of the steady-state 
output intensity lo on the pumping rate P^- To avoid erroneous contributions from high 
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FIG. 3: (Color online) Average steady-state emission intensity Iq (black x's) vs. pumping rate 
Pr for a random system with g = 1.0. Both Iq and Pr are plotted on logio scales to clearly 
show the regions of spontaneous emission (SE), amplified spontaneous emission (ASE), and laser 
emission (LE). Red dashed line and blue dotted line are linear fits to the intensities of SE (with 
nearly constant reabsorption at Pj. < 0.1) and LE (well above the lasing threshold at Pr > 10), 
respectively. The slopes, written next to the lines, are equal to one, refiecting linear increase of Iq 
with Pr- The intensity of ASE increases superlinearly with (-Pr)^- 

frequency components (mentioned in Sec. HT]), lo is found by a spectral integration 



where ki = ka — Aka = 27r/1.2 fim^^ and ku = ka + Aka = 2tt/0A ^m^^. 

Figure [3] plots log/o versus logP^. The total pumping rate is normalized such that at 
Pr = 1 (}ogiQ Pr = 0), the system without noise reaches the transparency point in the steady 
state {{psix))^ = 0). With noise, Pr = I is just below the transparency point {{paix))^ < 0), 
because the atomic population in level 2 (P22) is reduced by spontaneous emission. The 
spontaneous emission intensity is linearly proportional to p22- However, when P22 < Pu, the 
spontaneously emitted light can be reabsorbed. The amount of reabsorption is determined 
by Pii — P22, which varies with Pr. At very low pumping {Pr < 0.1), pn ^ P22, thus 
Pu ~ P22 — 1; and the amount of reabsorption is almost constant. As P22 increases linearly 
with Pr, the output intensity of spontaneous emission grows linearly with Pr. At higher 
pumping 0.1 < Pr < 1, the decrease of pn leads to a significant reduction in reabsorption. 
In fact the amount of reabsorption decreases nonlinearly with P^, resulting in a superlinear 
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FIG. 4: (Color online) Steady-state emission spectra with noise (upper black line in each 

panel) compared to those without noise (lower red line and crosses in the same panel) at the same 
pumping rate Pr for a random system with g = 1.0. The values of Pr are written in each panel. 



increase of Iq with _P,.. Once the pumping rate is large enough to induce a population 
inversion (pn < P22)i the spontaneously emitted light experiences a net amplification. The 
ASE intensity increases superlinearly with P^, as seen in Fig. |3l Even with the existence 
of population inversion, the rate of light amplification may be less than the leakage rate, 
and there is no lasing oscillation. Once the pump exceeds a threshold, light leakage is 
compensated by amplification, and lasing oscillation occurs. Well above the lasing threshold, 
the optical gain is saturated and the growth of Iq with becomes linear again. 

To ensure these results are not limited to the particular configuration considered here, 
the simulations are repeated with another random seed (to initialize the noise terms) and 
another realization of a random structure with the same g. The results are qualitatively 
similar. Slight differences arise due to stochasticity. 



B. ASE spectra 

We Fourier-transform the output fields to obtain the emission spectra. Figure H] shows 
the steady-state emission spectra with noise \E{k)\'^ in comparison to those without noise 
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for increasing pumping rates. At Pr = 1.00 [Fig. Hl^a)], there is no net gain. Without noise, 
the initial seed pulse dies away, and there is no signal at the steady state. With noise, the 
steady-state emission spectrum has a broad peak. It is centered at the atomic transition 
frequency ka = 10.5 fim~^, resembling the spontaneous emission spectrum. On top of it 
there are many fine spikes whose frequencies change chaotically from one time window of 
Fourier transform to the next. They result from the stochastic emission process with their 
spectral width determined by the temporal length of the Fourier transform. Above the 
transparency point at Pr = 1.02 [Fig. 111(b)], the broad emission peak grows and narrows 
spectrally. This behavior is typical of ASE. Since the optical gain is frequency dependent, 
the emission intensity closer to ka is amplified more than that away from ka, leading to 
a spectral narrowing. The stochastic emission spikes are also amplified, especially those 
closer to ka in frequency. As the pumping rate increases more [e.g., Pr = 1.04 in Fig. 
m^c)], the broad peak grows and narrows further. Without noise, the emission spectra are 
blank, since ASE is neglected. When Pr = 1.06 [Fig. 11(d)], a single peak appears in the 
emission spectrum without noise. This peak is a delta function with its "linewidth" merely 
determined by the integration time of the Fourier transformation. It shows lasing occurs in 
a single mode, which corresponds to the resonance of the passive system at = 11.6 fim~^ 
in Fig. [21 The lasing frequency is pulled towards ka at which gain is maximal. A further 
increase of Pr to 1.08 leads to lasing in a second mode that corresponds to the resonance 
at A; = 9.6 fim~^. Frequency pulling is also seen here. In the absence of noise, single mode 
lasing can be achieved by carefully adjusting Pr. This is no longer the case when noise is 
introduced. The emission spectra with noise look very different. There is clearly no single- 
mode lasing at any pumping rate. Intensity of the broad emission spectrum is modulated, 
as seen in Figs. Ill(d)-[ll(g). The emission intensities are enhanced not only at the frequencies 
of lasing peaks without noise, but also at some other frequencies. Optical amplification is 
stronger at the resonant frequencies of the system and narrows the resonance peaks that 
overlap spectrally without gain. Reduced overlap of resonance peaks with gain results in 
a spectral modulation of emission intensity. For Pr = 1.10 [Fig. 111(g)], it is possible to 
associate the three lasing peaks for the case without noise to resonance peaks with noise. 
Additional resonance peaks are also discernible for the case with noise. Unlike the stochastic 
emission spikes, the frequencies of resonance peaks are stable in time, although their heights 
may vary from one time window of Fourier transform to the next. Their spectral widths are 
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FIG. 5: (Color online) Spectral width 5k (black squares) and amplitude Ai (red circles) of the 
broad emission peak vs. pumping rate Pr- 

notably larger than those of the stochastic emission spikes. 

We extract the width of the broad emission spectrum in Fig. HJat different pumping rates. 
Because of the noisiness of the spectrum, we use a "Lorentz error function" to objectively 
obtain the spectral width. A Lorentzian function L{k) describes the spectrum 



where Ai is the amplitude, /cq is the center frequency, and 6k is the full-width at half- 
maximum (FWHM). The Lorentz error function, given by 



is used to fit the numerical data. The frequency range of the integration over k' is limited 
to ki = ka — Aka = 27r/1.2 fim~^ and ku = ka + Aka = 27r/0.4 fim~^. Because of the 
preferential amplification of light with frequencies closer to ka, the emission spectrum is 
narrowed around ka- Thus, the center frequency is fixed at ko = ka and not adjusted during 
the fitting js^] of Eq. (E]) to the data. Ai and 6k are the fitting parameters. 

The values of 6k and Ai obtained from the fitting is plotted against the pumping rate 
in Fig. O As Pr increases from 1.00 to 1.10, 6k first decreases rapidly then decreases 
more slowly. Ai displays a superlinear increase with Pr. Amplification allows the emission 
intensity to build up quickly around the atomic transition frequency. This increase results 
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FIG. 6: (Color online) (a) Steady-state emission spectra \E{k)\'^ with noise over a small frequency 
range illustrating stochastic emission spikes. Crosses mark the wavelengths at which the emission 
intensities are obtained by the Fourier transform. Arrows mark spikes identified by a three-point 
peak finding method, (b) Statistical distribution of frequency spacing of spikes P{5kn) at Pj- = 1.00 
(-|-'s) and Pr = 2.00 (x's), plotted on the linear scale (main panel) and logarithmic scale (inset). 
An exponential fit is marked in the inset by a straight red line. 



in a rapid narrowing of the emission spectrum. Such behavior has been seen experimentally 
{2!. Recent studies reveal that the spectral narrowing resembles a condensation process, 
as it can be predicted by a nonlinear differential equation identical to that governing the 
ultracold atoms 
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The stochastic emission spikes have also been observed in the ASE spectra experimentally 



[z-llOll. They are attributed to single spontaneous emission events which happen to take long 
open paths inside the amplifying random medium and pick up large gain. The emergence 
of these spikes does not rely on resonant feedback or coherent interference. Their spectral 
width is determined by the temporal duration of the emission pulse. In principle, our classical 
noise model does not account for spontaneous emission on the single photon level. However, 
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millions of photons are emitted and amplified in the macroscopic random media with gain. 
Thus the quantum nature of photons can be ignored. We found the stochastic spikes in the 
emission spectra of our simulation bear similar characteristics to the ASE spikes measured 
experimentally. An example of the stochastic spikes is shown in Fig. M^a), which is an 
enlargement of the emission spectrum in Fig. ID^a). We find that the spectral width of the 
stochastic spikes is determined by the integration time of the Fourier transformation Tp. 
Because of the long Tp, the spikes are usually much narrower than the lasing peaks, as long 
as the pumping rate is not too high. In the emission spectrum the intensity is calculated at 
the frequency step 6kn that is determined by Tp. If the intensity at k is larger than those at 
k ± 6kn, a spike is found at k. Using this three-point peak-finding method, we extract the 
frequencies of spikes from the calculated spectra as shown in Fig. M^a), and compute the 
frequency spacing of adjacent spikes. Figure [6](b) plots the statistical distribution P{6kn) of 
frequency spacing between adjacent stochastic spikes for = 1.00 and Pr = 2.00. The two 
distributions coincide, revealing P{6kn) is independent of Pr- As evident from the log-linear 
plot in the inset of Fig. [6](b), P{6kn) decays exponentially at large 6kn- This behavior is 
identical in character to the experimental result [l^. Note that the level-off of P{6kn) at 
smaller 6kn is an artifact of limited spectral resolution. 

C. Lasing modes 

The resonance peaks, which are hardly seen in the emission spectra for Pr < 1.10 (Fig. 
4), grow rapidly as Pr increases further above 1.10 (Fig. [7]). They become narrower and 
well separated, surpassing the stochastic emission spikes. There are clearly more peaks 
in the emission spectra with noise than those without noise at the same pumping level. 
This is because all modes are constantly excited by the noise and subsequently amplified 
in the presence of population inversion. Hence, the pump energy is distributed over more 
peaks. Nevertheless, all the lasing peaks without noise correspond to strong emission peaks 
with noise. We enumerate six major peaks in Fig. El^d), with 1 being the strongest. For 
Pr > 2.00 [Figs. Uli) and[71^g)], the difference between the emission spectra with noise and 
those without noise is reduced. 

We compare the intensities and frequencies of the lasing modes with noise to those without 
noise. In both cases, the strongest peaks are 1 and 2 (Fig. [7j). We plot their intensities 
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FIG. 7: (Color online) Steady-state emission spectra with noise (upper black line in each 

panel) compared to those without noise (lower red line and crosses in the same panel) at the same 
pumping rate Pr for a random system with g = 1.0. The values of Pr are written in each panel. 

versus Pr in Fig. [8]^a). In the absence of noise, there is a clear threshold for lasing. For 
example, mode 1 has zero intensity for P,. < 1.06. Once Pr exceeds 1.06, its intensity rises 
quickly. The sharp turn-on at Pr = 1.06 marks the lasing threshold for mode 1. Mode 2 
reaches its lasing threshold by Pj. = 1.08 and its intensity increases almost linearly with Pr. 
Although modes 1 and 2 display notable frequency pulling just above the lasing threshold, 
their frequencies do not shift significantly as Pr increases further above the threshold. This 
is due to gain saturation. At Pr > 1.10, the center frequencies of lasing peaks with noise 
are almost identical to those without noise. Thus noise does not affect the frequencies of 
lasing modes. Since each lasing peak has a finite width, we integrate the emission intensity 
over a spectral range set by the mid frequencies between adjacent peaks. As shown in 
Fig. El^a), the intensities of modes 1 and 2 increase gradually with Pr. The soft turn-on 
makes it difficult to pinpoint the exact value of the lasing thresholds. Because of ASE, the 
modal intensity is non-zero below the lasing threshold. Above the threshold pumping rate 
for lasing without noise, the intensity with noise is notably lower than that without noise 
because some pump energy is diverted to other modes via ASE. The superlinear increase 
of modal intensity around the threshold is caused by ASE. It is evident that noise almost 
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FIG. 8: (Color online) Modal intensities with noise (solid lines) compared to those without noise 
(dashed lines) for a random system with g = 1.0. (a) Intensity of mode 1 Ii (red thick lines) and 
mode 2 I2 (blue thin lines) vs. Pr- (b) Intensity of mode 3 with (green thick lines) and mode 4 
/4 (black thin lines) vs. Pr- Modes are enumerated in Fig. [Tl^d). 



equalizes the intensities of modes 1 and 2, despite the fact that mode 1 is clearly stronger 
than mode 2 without noise. 

Figure Et^b) shows the intensity of modes 3 and 4 [enumerated in Fig. El^d)] with and 
without noise as increases. Without noise, mode 3 reaches its lasing threshold at = 
1.08. Mode 4 has a similar lasing threshold, but its amplitude remains small until Pr = 1.22 
[Fig. Cl^b)]. With noise, the intensities of both modes start rising from zero at Pr < 1.08. 
They increase superlinearly with Pr and are greater than the intensities without noise even 
about the threshold for a small range of pumping rates [inset of Fig. El^b)]. This is most 
noticeable for mode 4 in the range Pr < 1.26. The co-existence of multiple modes and 
their interactions through the gain material make it difficult to define the lasing threshold 
for each mode using previously developed methods for single mode lasers 40l-l44|. Though 
the lasing threshold is not precisely defined here, the soft turn-on and subsequently smaller 
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an -1 
o 



slope = -0.97 +/- 4.6% 
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FIG. 9: (Color online) Linewidths of lasing mode 1 5ki (red crosses) and lasing mode 2 Sk2 (blue 
diamonds) versus the corresponding mode intensity Ii and I2 for a random system with g = 1.0. 
Linear fits to the data give the power by which the linewidths decrease. Mode 1 oc I^^'^"^ . Mode 2 
oc 

intensities at larger Pr in the case with noise shows the threshold is increased for each of 
the four dominant modes (1-4) when noise is included. 

Next we calculate the spectral width of lasing modes, that is impossible to do with the 
noiseless simulation. Considering the noisiness of the spectrum, we again use a Lorentz error 
function [Eq. (^] to obtain the linewidth objectively. The integration of emission intensity 
for the Lorentz error function is limited to the spectral range of each mode, which is the 
same as that used to obtain the spectrally-integrated intensity. Figure [9] plots the linewidths 
6k of modes 1 and 2 with respect to the steady-state intensities /. Mode 1 narrows the most 
dramatically; its linewidth decays over two orders of magnitude. On a log-log scale, the 
data for each mode falls onto a straight line, indicating a power-law decay. We fit the data 
by 6k oc /° within a range Ii < I < lu- /; is set by the threshold pumping rate without 
noise, at which separate resonance peaks emerge in the presence of noise. /„ is determined 
by the pumping rate at which an accurate estimate of the linewidth is no longer possible 
due to limited spectral resolution (determined by the running time of the simulation). For 
mode 1, the exponent a = —0.97 ±4.6% is close to the Schawlow-Townes prediction of laser 



linewidth 



45| . For mode 2, a = —0.70 ±4.4% so the linewidth decays slower, probably due 



;o mode competition for gain. It is known that multimode operation affects laser linewidths 
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391]. A more quantitative investigation will be carried out in the future. 



Finally we look at some of the smaller peaks, e.g., 5 and 6 in Fig. El^d). Both have 
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corresponding peaks in the noiseless spectra, but they are orders of magnitude smaller than 
the main peaks and cannot be seen on the vertical scale of Fig. [71 With noise, peak 6 
has a much larger amplitude and is visible together with the major peaks in the emission 
spectrum. Close by the frequency of peak 5, there are two resonances, one on either side of 
ka in Fig. |2]^b). At lower pumping rates, only a "composite" peak appears at k = 10.3 /im~^. 
The linewidths of the two modes exceed their frequency spacing, which is reduced by the 
frequency pulling effect. Consequently, the two modes are indistinguishable and appear to 
be merged. At higher pumping rates, their linewidths decrease further, but the amphtudes 
remain relatively small compared to the four main peaks. Such weaker modes are affected 
more by gain nonlinearity, and display complicated behavior with increasing pumping rate. 
Detailed investigation of this behavior will be left for future studies. 

VI. ASE AND LASING IN A RANDOM SYSTEM WITH NON-OVERLAPPING 
RESONANCES 

In this section we study laser emission characteristics of the ID random system with 
g = 0.5. With higher refractive index contrast (An = 0.25), light leakage is reduced and 
so is the lasing threshold. Figure shows the steady-state emission spectra for increasing 
pumping rates with and without noise. Pr is normalized to the value at which = in the 
absence of noise. In Fig. [TUT a). Pr = 1 and there is no gain, so without noise the steady-state 
emission intensity is zero. With noise, the steady-state emission spectrum has a broad peak 
around the atomic transition frequency. Spectral modulation of emission intensity is evident. 
Higher emission intensities match the transmission peaks in Fig. [T| lower intensities match 
the transmission dips. Because the Thouless number is less than unity, the quasimodes are 
already separated. Particular modes may be even narrower and farther apart. They appear 
as peaks in the emission spectrum even without gain [Pj. = 1.0). 

With only a slight increase of the pumping rate to Pr = 1.02 [Fig. [TOT b)]. a single lasing 
peak appears at = 10.8 fim"^ in the absence of noise. With noise present, the broad 
emission peak grows and narrows around ka- Intensity modulation is enhanced, as the 
resonance peaks become narrower by light amplification. For P,. > 1.04 [Figs. fTor c)- [Tor g)]. 
well separated peaks develop in the emission spectra with noise. The two strongest emission 
peaks, enumerated in Fig. [TOlfd). have the same frequencies as the lasing peaks without 
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FIG. 10: (Color online) Steady-state emission spectra with noise (upper black line in each 

panel) compared to those without noise (lower red line and crosses in the same panel) at the same 
pumping rate Pr for a random system with g = 0.5. The values of Pr are written in each panel. 



noise. They correspond well to the two resonances nearest ka in the passive system (boxed 
in Fig. |2]). By = 2.00 [Fig. [TOT g)] the spectrum with noise resembles that without noise. 
In both cases, the number of major peaks is three. The influence of noise is reduced for the 
system of smaller g, because the lower lasing threshold narrows the range of pumping rates 
where ASE dominates. With increasing Pr, gain saturation quickly sets in to suppress the 
fluctuations. 

Figure [11] plots the spectrally-integrated intensity of the two peaks enumerated in Fig. 
[TO](d). Compared to the two modes in Fig. [8]^a), the increase of intensity with pumping is 
more rapid. As before, we cannot pinpoint the exact lasing threshold for each mode because 
of multimode operation. Nevertheless, it is evident that without noise the onset of lasing 
oscillation occurs at a lower pumping rate, and the modal intensity is higher than that 
with noise. This is because the pump energy is partly consumed by ASE in other modes 
in the presence of noise. However, by Pr = 2.00 the effect of noise is diminishing, and the 
modal intensity with noise converges to that without noise. Hence, the transition from the 
ampliflcation of spontaneous emission to the lasing oscillation happens over a relatively small 
range of pumping rates, and the effect of noise is less signiflcant than that in the random 
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FIG. 11: (Color online) Modal intensities with noise (solid lines) compared to those without noise 
(dashed lines). Intensity of mode 1 (red thick lines) and mode 2 (blue thin lines) vs. Pr for a 
random system with g = 0.5. Modes are enumerated in Fig. [TUTd). 



system with larger g. 



VII. CONCLUSION 



The effects of fluctuations caused by interactions of atoms with reservoirs were studied in 
random lasers. A FDTD-based method for solving the stochastic Maxwell-Bloch equations 
was employed. It is particularly well-suited for studies of light-matter interaction in complex 
systems without prior knowledge of resonant modes. Two random systems with different 
degrees of spectral overlap of resonances were investigated. We were able to simulate ampli- 
fied spontaneous emission (ASE) below the lasing threshold and capture the transition from 
ASE to lasing. 

In the case of overlapping resonances, the emission spectra at low pumping are broad 
and continuous. Above the transparency point, frequency-selective amplification leads to a 
dramatic narrowing of the emission spectrum and a superlinear increase of the peak emission 
intensity. Such behavior is in accordance with early experimental results {2]. Moreover, our 
simulation reproduced the stochastic emission spikes in the spectra, with similar character- 

nn 

istics to the experimentally observed ASE spikes [9|, 110|. Previous experiments found the 
spectral width of ASE spikes depends on the temporal duration of the emission pulse. Here, 
we found the width of stochastic emission spikes is determined by the integration time of 
the Fourier transform of the output field. The statistical distribution of frequency spacing 
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of spikes displays an exponential tail, as seen experimentally. The spikes have no relation to 
the resonant modes of the system, and can be clearly differentiated from the emission peaks 
formed by resonances. 

We compared the lasing behavior with noise to that without noise in the same random 
system. The lasing peaks in the spectra without noise coincide with peaks in the spectra with 
noise. Hence, noise does not affect mode frequencies. However, all modes within the gain 
curve are constantly excited by noise and subsequently amplified by stimulated emission. 
Therefore, there are always multiple modes appearing in the steady-state emission spectra. 
The regime of single mode lasing, realized in the noiseless simulation by fine tuning of pump, 
disappears. With some portion of pump energy diverted to ASE in other modes, the lasing 
modes have higher thresholds than those without noise. Moreover, the ASE below the lasing 
threshold results in a soft turn-on of the lasing mode. It is in sharp contrast to the abrupt 
turn-on in the simulation without noise, where the emission intensity vanishes below the 
lasing threshold. When the pumping rate is well above the threshold value, the spectra 
with noise become more similar to the spectra without noise, both showing multimode 
lasing. Thus, noise has the greatest influence on lasing behavior near threshold. With noise 
included, we can calculate the spectral widths of individual lasing modes, and observe their 
decrease with increasing pump. The decrease appears to follow the Schawlow-Townes law 
for the strongest lasing mode, but not for other modes, probably due to mode interactions. 

The effects of noise on lasing become less significant in random systems with a smaller 
degree of spectral overlap of resonances. If the Thouless number is less than unity, the 
resonant modes can be resolved in the emission spectra below the transparency point. ASE 
narrows the resonance peaks, making them more distinct. The transition from ASE to lasing 
occurs over a narrower range of pumping rate, because the lasing threshold is lower in the 
case of non-overlapping modes. The difference between the simulation results with less mode 
overlap and those with larger mode overlap agrees qualitatively to the experimental data 
23| that compare different particle densities. Increasing the refractive index contrast An in 
our simulation enhances the scattering strength, which is similar to increasing the density 
of scattering particles in the experiments. 

These studies shed light on the transition from ASE to lasing in random systems, that 
is poorly understood. The results presented here are limited to the steady state. Noise is 
expected to have a greater effect on the dynamics, e.g., the buildup of lasing modes and 
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temporal fluctuations and switching of lasing modes on short time scales. These phenomena 
can be studied with our numerical method. Furthermore, this method can be extended to 
the study of random lasing in higher dimensions. The larger density of modes and potentially 

stronger mode overlap in frequency may enhance the noise effects. 
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